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We use the recently proposed In-Medium Similarity Renormalization Group (IM-SRG) to carry 
out a systematic study of closed-shell nuclei up to 56 Ni, based on chiral two- plus three- nucleon 
interactions. We analyze the capabilities of the IM-SRG by comparing our results for the ground- 
state energy to Coupled Cluster calculations, as well as to quasi-exact results from the Importance- 
Truncated No-Core Shell Model. Using chiral two- plus three-nucleon Hamiltonians whose resolution 
scales are lowered by free-space SRG evolution, we obtain good agreement with experimental binding 
energies in 4 He and the closed-shell oxygen isotopes, while the calcium and nickel isotopes are 
somewhat overbound. 
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I. INTRODUCTION 

The adaptation of Renormalization Group (RG) and 
Effective Field Theory (EFT) methods has led to many 
advances in the treatment of the nuclear many-body 
problem in recent years 0, 0] . The Similarity Renormal- 
ization Group (SRG), in particular, has become the tool 
of choice to lower the resolution scale of two- (NN) and 
three-nucleon (3 N ) interactions [U-IH • The improved con- 
vergence properties of the resulting effective interactions 
make them an ideal input for many ab initio methods. As 
a consequence, the range of nuclei to which such methods 
can be applied successfully has increased significantly. 

Recently, the In-Medium SRG (IM-SRG) has been pro- 
posed as a new ab initio many-body technique, which 
uses the powerful SRG flow-equation framework to di- 
rectly calculate nuclear structure observables, or to de- 
rive effective Hamiltonians for Shell Model calculations, 
which can provide complete spectroscopic information for 
closed- and open-shell nuclei from first principles @, 0] ■ 
Here, we use the natural capabilities of the IM-SRG to 
include 3iV interactions in normal-ordered two-body ap- 
proximation, and present a systematic study of closed- 
shell nuclei up to 56 Ni, based on chiral NN and 37V in- 
teractions. 

In Sect. HU we give an overview of the IM-SRG formal- 
ism, and discuss its connection to Many-Body Perturba- 
tion Theory (MBPT). Section [m] describes the Hamil- 
tonians used in this work and provides some technical 
details on our calculations. In Sect. IIVI we compare IM- 
SRG ground-state energies with results from the non- 
perturbative Coupled Cluster (CC) method at different 
truncation levels (see, e.g. §), as well as the Importance- 
Truncated No-Core Shell Model (IT-NCSM) (UGH. The 
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IT-NCSM approach is a diagonalization of the Hamilto- 
nian in a Hilbcrt space spanned by importance-sampled 
basis states. Aside from small uncertainties associated 
with the importance sampling and model space extrapo- 
lations, the IT-NCSM results are exact, and particularly 
useful for quantifying the uncertainties of our IM-SRG 
results for closed-shell nuclei, which are presented at the 
end of Sect. HVl 



II. FORMALISM 

A. Normal Ordering 

Our starting point is an intrinsic nuclear A-body 
Hamiltonian containing both NN and 3iV interactions, 



H 



(1) 



where is the two-body part of the intrinsic kinetic 
energy (see, e.g., [IH). Choosing a single Slater deter- 
minant |(f>) as the reference state, we can rewrite the 
Hamiltonian in terms of normal-ordered operators, 
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where the indices include all single-particle quantum 
numbers, and the strings of creation and annihilation op- 
erators obey 
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It is convenient to work in the eigenbasis of the one-body 
density matrix in the following, so that 



(= in n 



2 



and the individual normal-ordered contributions in 
Eq. (3J) are 
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Due to the occupation numbers in Eqs. ©-(CO), the sums 
run over occupied (hole) states only. Note that the zero-, 
one-, and two-body parts of the Hamiltonian all contain 
in-medium contributions from the free-space 37V inter- 
action. The normal-ordered 37V contribution W will be 
omitted in the following, leading to the normal-ordered 
two-body approximation (N02B), which has been shown 
to be a very good ap proximation for the nuclei considered 
in this work 



B. IM-SRG Flow Equations 

The aim of the IM-SRG is to decouple the ground- 
state of the Hamiltonian from all excitations by means of 
a continuous unitary transformation. The transformed 
Hamiltonian is defined as 



H(s) = U J '{s)H(0)U{s), 



(9) 



which, upon taking the derivative w.r.t. the flow parame- 
ter s, yields the following first-order operator differential 
equation: 



-H(s) = [r)(s),H(s)], 
as 

with the generator formally defined by 
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When carried out exactly, the IM-SRG is a unitary 
transformation in A-nucleon space, and consequently, 
r/(s) and H(s) are A-body operators. When Eq. (JTUJ) is 
integrated, every evaluation of the commutator increases 
the particle rank of H(s), e.g., 
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All of these induced contributions will in turn contribute 
to the parts of H(s) with lower particle rank in subse- 
quent integration steps. Because an explicit treatment of 
all contributions up to the ^4-body level is clearly not fea- 
sible, we have to introduce a truncation to close the sys- 
tem of IM-SRG flow equations. The simplest approach 
is to truncate H(s) at a given particle rank n < A, 
which is motivated by the cluster decomposition prin- 
ciple for short-range interactions (see, e.g., [3]). It has 
been shown that the omission of W, the residual 3iV part 
of the Hamiltonian (|2|), is a good approximation as long 
as the in-medium contributions of the free-space 3iV in- 
teraction are accounted for by the normal-ordered zero-, 
one-, and two-body interactions [l^ - fl^ |. We therefore 
truncate both H(s) and rj(s) at the two-body level, and 
refer to this truncation as IM-SRG(2). An additional 
truncation which was motivated by a perturbative anal- 
ysis was proposed in Q , but will not be considered in this 
work. For the Hamiltonians considered in the following, 
the numerical results of the two truncation schemes agree 
to within a few keV @. 

Note that the presence of the commutator in the 
flow equation (fT0)l guarantees the IM-SRG to be size- 
extensive, i.e., the IM-SRG wave-function U(s) 1$) can 
be expanded in terms of linked diagrams only [R5H18| . 
As a consequence, errors which are introduced by trun- 
cating the many-body expansion scale linearly with A 
[ill [l8|| , and we can, in principle, judge the quality of the 
IM-SRG results in medium- or heavy-mass nuclei based 
on comparisons with results from exact methods like the 
No-Core Shell Model (NCSM) [TJ, S3, which are only 
available for light nuclei. 

Since the focus of this work is on closed-shell nuclei, we 
assume spherical symmetry. Single-particle indices col- 
lectively represent the radial, angular momentum, and 
isospin quantum numbers i = (kilijiTi), and do not de- 
pend on the angular momentum projection m,. The ma- 
trix elements of single-particle operators are diagonal in 
all but the radial quantum numbers, e.g., 
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and two-body matrix elements are coupled to good J 
and independent of M . Suppressing all s-dependence for 
brevity, the resulting J-scheme flow equations read: 
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where j = \/2j + 1, indices with a bar indicate time- 
reversed states, and the rj and T matrix elements in the 
last line of Eq. ([To]) are obtained by a generalized Pandya 
transform (see, e.g., plj). 
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C. Choice of Generator 

The next step is to identify the "off-diagonal" part of 
the Hamiltonian, which we wish to suppress by integrat- 
ing the flow equation. To this end, we consider the repre- 
sentation of the flowing Hamiltonian H(s), truncated to 
two-body operators, in a many-body basis consisting of 
up to A-particle-A-hole (ApAh) excitations of the refer- 
ence state |<I>), as sketched in Fig.[TJ The OpOh reference 
state is coupled to lplh and 2p2h excitations by the ma- 



OpOh lplh 2p2h 3p3h 



OpOh lplh 2p2h 3p3h 



(i\ H(0) \j) 



(i\H(oo)\j) 



FIG. 1. Schematic representation of the initial and final 
Hamiltonians, H(0) and H(oo), in the many-body Hilbert 
space spanned by particle-hole excitations of the reference 
state. 
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and their Hermitian conjugates, which define the one- 
and two-body pieces of the off-diagonal Hamiltonian H od . 
This H od is now used to construct the generator 0, [22j , 
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The energy denominators appearing in Eq. (|20|) represent 
a link between the IM-SRG and Many-Body Perturba- 
tion Theory (MBPT), which will be discussed in Section 
III Dl As suggested by White in Ref. [22|, we use the 
Epstcin-Nesbet partitioning rather than the more com- 
mon M0ller-Plesset partitioning [H and construct energy 
denominators using the diagonal matrix elements of H in 
our chosen matrix representation (Fig.[T]). This naturally 
regularizes the generator in situations where the differ- 
ence of the single-particle energies might become small 
during the IM-SRG flow. 



*-Q*<Q> O-<0>-<8>*<8>*<8>*- 

FIG. 2. Schematic illustration of the energy flow equation (|14|) for the White generator (|20[) in terms of Hugenholtz diagrams 
(see text). The grey vertices correspond to H(s), and the double lines indicate energy denominators calculated with f(s). On 
the right-hand side, the flow equation is expanded in terms of H(s — 5s) (simple black vertex) and the corresponding energy 
denominators from f(s — 5s) (single lines). The dots indicate omitted terms from the expansion of the Epstein-Nesbet energy 
denominator, as well as additional diagrams generated by the integration step s — 5s — > s. 



Because rj is now given by ratios of energies, / and T 
only contribute linearly to the magnitude of the right- 
hand side of the IM-SRG flow equations ([I4] l -([T6 |) . This 
reduces the stiffness of the flow equations significantly 
compared to calculations with the canonical Wegner gen- 
erator, rf = [H, H od ] HHH, where third powers of / and 
r appear. 

The use of (f20|) also implies that the off-diagonal ma- 
trix elements (fT5|) , (fTO|) are suppressed at all energy scales 
simultaneously, i.e., 

H od (s) ~ H od (0)e- s , (22) 

in contrast to the canonical Wegner generator, which be- 
haves in a true RG sense and preferentially suppresses 
matrix elements between states with large energy differ- 
ences (see the discussion for the free-space SRG in [l[). 
However, since our aim is to solve the many-body prob- 
lem via the IM-SRG transformation, we are only inter- 
ested in the limit s — > 00, and then the White generator 
(|20| must yield the same result as the Wegner generator, 
up to truncation errors. We have verified numerically 
that the energies obtained with these generators differ at 
the single- keV level in the present work (also see Ref. 

!)• 



of the second-order MBPT correction to the ground- 
state energy, but note that / and T both depend on the 
flow parameter s. Thus, the ground-state energy is RG- 
improved by including contributions from higher-order 
diagrams into E(s). 

The effect of integrating the energy flow equation (fT4|) 
by a single step s — 5 s — > s is illustrated schemati- 
cally in Fig. [2] Second-order terms in H(s), which are 
added to E(s), contain contributions from second- and 
third-order diagrams when expressed in H(s — 5s). Note 
that all topologies for third-order diagrams are gener- 
ated: the second line of the two-body flow equation (fTB)) 
adds particle-particle and hole-hole ladder contributions 
of r(s — 5s) to r(s), while the third line adds particle- 
hole contributions, and upon insertion of T(s) into the en- 
ergy flow equation, we obtain the corresponding Hugen- 
holtz energy diagrams (see, e.g., [1, HH) in Fig. [2] Thus, 
integrating the IM-SRG flow equations from s = to 
00 corresponds to performing a partial re-summation of 
the MBPT series. This includes energy corrections for 
the initial normal-ordered Hamiltonian ([2]) which are (at 
least) complete to third order, as well as ladder, ring, and 
classes of ladder-ring-interference diagrams to all orders. 
A detailed diagrammatic analysis of the IM-SRG will be 
presented in a forthcoming publication [25| . 



D. Connection to Perturbation Theory 



III. HAMILTONIANS AND IMPLEMENTATION 



In the previous section, we have already indicated that 
the White generator (|2TH) reveals a clear connection be- 
tween IM-SRG and MBPT. If we substitute Eq. fl2DJ) into 
the flow equation for the ground-state energy, Eq. (fT4|) . 
we obtain 
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where we have expanded the Epstein-Nesbet energy de- 
nominator as a Neumann series to make the connection 
with the more commonly used M0ller-Plesset partition- 
ing clearer Q. Equation (f2"3"|) clearly has the structure 



In this work, we consider two classes of Hamiltoni- 
ans. The first one, referred to as NN + 37V-induced, 
is based on a chiral NN interaction at next-to-next-to- 
next- to- leading order (N 3 LO) by Entem and Machlcidt, 
with cutoff A = 500 MeV [26|. The resolution scale of 
this Hamiltonian is subsequently lowered by means of a 
free-space SRG evolution 0, H3, EH , and the induced 3JV 
interactions are kept explicitly. 

The second class of Hamiltonians is referred to as 
A^A^+SiV-full in the following, and supplements the afore- 
mentioned initial NN Hamiltonian with a local chiral 3iV 
interaction at next-to-next-to-leading order (N 2 LO) [30| . 
For a cutoff A = 500 MeV, this Hamiltoni an p erforms 
well in calculations for s- and p-shcll nuclei [20l . \2§L [3l| , 
but when it is evolved to lower resolution scale in order 
to improve convergence for calculations in the upper p- 
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FIG. 3. (Color online) Convergence of the IM-SRG(2) ground-state energy of closed-shell nuclei for the NN + 3iV-full Hamil- 
tonian (see text), evolved to A = 2.0 fm _1 . The energies indicated by dashed lines are obtained with the extrapolation method 
proposed in [2S|. 



shell and beyond, substantial 47V interactions are induced 
[H HI]. Following Ref. [H, we reduce the cutoff of the 
initial 3N force to A = 400 MeV, and use the low-energy 
constants co = —0.2 and ce = 0.098 that have been 
refitted to the He binding energy and triton beta-decay. 

The SRG evolution of both types of Hamiltonians is 
carried out in two- and three-body Jacobi harmonic oscil- 
lator (HO) bases with large but finite dimension [27l |32|| . 
One has to ensure that this model space captures all 
relevant contributions of the initial Hamiltonian. For 
small HO frequencies, in particular, the truncation of the 
model space at fixed HO excitation energy, up to AQMl 
in the present study, leads to additional uncertainties 
(see [32] for a detailed discussion). For this reason, the 
smallest oscillator parameters we use in the following are 
Ml = 20 MeV in the oxygen isotopes, and Ml = 24 MeV 
in heavier nuclei. This can be avoided in future appli- 
cations by either improving the SRG model space in the 
HO basis [32| or by carrying out the 3N evolution in 
momentum space [5j. 

The SRG-evolved NN + 3N Hamiltonian is subse- 
quently transformed to an angular-momentum coupled 
basis built from single-particle spherical HO states with 
quantum numbers e = 2n + l < eMax- For the 3iV Hamil- 
tonian, we introduce an additional cut e\ + 62 + e% < 



< 3eMax to manage the storage requirements of 
the matrix elements. Throughout this work, we use 
£^3Max = 14. As shown in Ref. [ljj, p/-shell nuclei 
are currently the upper limit of ab initio calculations 
with 3A^ interactions, because interactions between nu- 
clcons in higher shells are strongly affected or even en- 
tirely removed by the i?3Max cut. In actual IM-SRG(2) 
calculations, we find that increasing i?3Max from 10 to 
12 and from 12 to 14 changes the ground-state ener- 
gies by 2% and 1%, respectively, in the calcium iso- 
topes, and less so in lighter nuclei, for resolution scales 
A = 1.88, . . . , 2.24 fm _1 . For larger A, the resulting un- 
certainty of the ground-state energies grows significantly. 

To obtain reference states, we solve the Hartree-Fock 
equations using the code described in [33j], which has 
been extended to include 3A^ forces. The full three-body 
Hamiltonian is used in our calculations at this stage. 
Once a converged HF ground-state is obtained, the nu- 
clear Hamiltonian is normal-ordered w.r.t. to this solu- 
tion, and the resulting in- medium zero-, one-, and two- 
body operators serve as the initial values for the IM-SRG 
flow equations. As discussed above, the residual three- 
body term is neglected. The flow equations are subse- 
quently integrated using the solver CVODE [34| until 
decoupling is achieved. To check the convergence, we 
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calculate the energy correction from second-order MBPT 
for the flowing Hamiltonian H(s), which is entirely due 
to the off-diagonal part of the Hamiltonian as defined in 
Eqs. (HHJ) and (fTO)l . We assume that sufficient decoupling 
is achieved once the perturbative contribution is smaller 
than 10~ 6 McV, corresponding to relative changes in the 
flowing ground-state energy of 10~ 7 or less in medium- 
mass nuclei. 

To illustrate the convergence of our calculations, we 
show IM-SRG(2) ground-state energies as a function of 
the model space size eMax and the oscillator parame- 
ter Ml for the SRG-evolved TVTV + 3TV-full Hamiltonian 
at A = 2.0 fm _1 in Fig. [3] The rapid convergence of 
these ground-state energies is also representative for the 
TVTV + 3V-induced Hamiltonian in the range of studied 
resolution scales A = 1.88, 2.24 fm _1 . The dashed 
lines in Fig. [3] indicate energies that were obtained with 
an extrapolation to infinite HO spaces as proposed in 
[29j . Note that the p/-shell nuclei, in particular 56 Ni, 
exhibit kinks in the energies around Ml = 24 MeV which 
are caused by the free-space SRG model space trunca- 
tion discussed above. We will come back to this subject 
when we summarize our ground-state energies at the end 
of Sect. ED 



IV. RESULTS 

In this section, we discuss IM-SRG results for the 
ground-state energies of closed-shell nuclei. We first show 
how the IM-SRG re-sums correlations into the ground- 
state energy as we integrate the flow equations. Next, 
we compare IM-SRG results to other non-perturbative 
methods to assess its capabilities. Trends of the ground- 
state energies under the variation of the resolution scale 
A help us to identify and quantify the uncertainties of our 
results. We conclude the discussion by giving an overview 
of IM-SRG ground-state energies for closed-shell nuclei 
up to 56 Ni. 

Let us now demonstrate how the IM-SRG re-sums 
the MBPT series, as explained in Section III Dl To this 
end, we show the flowing 40 Ca ground-state energy plus 
second- and third-order MBPT corrections as a function 
of the flow parameter s in Fig.[4j using both the bare NN 
Hamiltonian and the TV TV + 3TV-induced Hamiltonian at 
a resolution scale A = 2.0 fm _1 . Evidently, correlations 
are summed very rapidly into E(s), and after only 25-30 
integration steps (s ~ 2 — 3), the size of MBPT correc- 
tions is reduced to less than 1% of E(s). Note, however, 
that the initial size as well as the efficiency of the re- 
summation strongly depends on the initial Hamiltonian. 
For the bare TVTV Hamiltonian (see Sect. IIII|) . the refer- 
ence state is unbound at E(0) = 114 MeV and only the 
expansion of single-particle wave-functions in terms of 
HO states keeps the nuclei from breaking apart. The ini- 
tial second- and third-order contribution are —434.1 MeV 
and —14.2 MeV, respectively, but by comparing with 
the final IM-SRG result E = -349.6 MeV. we find 
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FIG. 4. (Color online) IM-SRG (2) ground-state energy of 
40 Ca as a function of the flow parameter s for the bare TVTV 
Hamiltonian (top) and the TVTV + 3TV-induced Hamiltonian 
with A = 2.0 fm" 1 (bottom). Shown are E(s) (•), and E(s) 
plus second (■) and third-order (♦) MBPT corrections for 
H(s). The dashed line indicates the final value E(oo). All 
calculations were done with eMax = 10, at optimal Ml = 
32 MeV (top) and 24 MeV (bottom). For s > 1.5 we only 
show every tenth numerical data point to reduce clutter. 



that the re-summed higher-order contributions are about 
— 16 MeV, i.e., larger than the third-order term. This in- 
dicates the poor convergence properties of the MBPT 
series. For the soft TVTV + 3TV-induced Hamiltonian with 
A = 2.0 fin -1 (see Sect. ILU) . on the other hand, HF + 
second- and third-order MBPT is a good approximation 
to the ground-state energy, and the higher-order correc- 
tions included via the IM-SRG only amount to about 
—500 keV. Note, however, that even for soft Hamilto- 
nians, partial summations of the MBPT series to finite 
order do not converge systematically, as demonstrated in 

[nm. 

Next, we compare our IM-SRG(2) results for the 
ground-state energies of 4 He, 16 0, and 40 Ca with other 
non-perturbative methods. To this end, we consider 
energies from Coupled Cluster calculations with single 
and double excitations (CCSD) and A-CCSD(T) to as- 
sess contributions from triple excitations @, EH EH ■ The 
energies of both methods are obtained using the N02B 
approximation (see Sect. |Hj and i?3Max = 14, i.e., the 
CC and IM-SRG results are based on exactly the same 
input matrix elements. As a benchmark, we use ener- 
gies from the quasi-exact Im por tance- Truncated No-Core 
Shell Model (IT-NCSM) 0,ll3|, which contain the com- 
plete 3V interaction. We note that the i?3Max cut is 
naturally compatible with the model-space truncation of 
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FIG. 5. (Color online) Ground-state energies of stable closed- 
shell nuclei as a function of the resolution scale A for the 
SRG-evolved NN + 3AT-induced (left column) and NN + 3N- 
full Hamiltonians (right column), using different many-body 
methods: IM-SRG(2) (•), CCSD (■), A-CCSD(T) (♦), and 
IT-NCSM (a). IM-SRG and CC results are obtained with 
eMax = 14, and the IT-NCSM results have been extrapolated 
to infinite model space, with error bars indicating the uncer- 
tainties of the importance truncation and the model space 
extrapolation (for 4 He, the error bars are smaller than the 
symbols). The gray solid line for 4 He in the top left panel 
is the result of a converged NCSM calculation with the bare 
NN Hamiltonian, E = -25.39 MeV (see e.g. Q). Dashed 
black lines are experimental values from [37j . 



the IT-NCSM, and therefore not a source of additional 
uncertainties in IT-NCSM calculations. In Fig. [5j we 
show such a comparison for the NN + 3-/V-induced and 
NN + 37V-full Hamiltonians at various resolution scales 
A. 

Surveying the results, we observe that as A decreases, 
the differences between the IM-SRG(2), CCSD and A- 
CCSD(T) energies becomes smaller, and all three meth- 
ods approach the quasi-exact IT-NCSM results in 4 He 
and 16 O. The only deviation from a fairly uniform pat- 
tern occurs for 4 He with the NN + 3A^-full Hamiltonian, 
where the difference between CCSD and IT-NCSM seems 
to increase slightly, but here the absolute scales of the 
energy are quite small. The improving agreement of IM- 





-320- 



-340- 



-360- 



FIG. 6. (Color online) Ground-state energies of stable closed- 
shell nuclei for the bare NN Hamiltonian (A = oo), using IM- 
SRG(2) (•), CCSD (■), and A-CCSD(T) (♦). IM-SRG and 
CC results are obtained with eMax = 14. The gray solid line 
for 4 He indicates the result of a converged NCSM calculation, 
E = -25.39 MeV (see e.g. H). The open and solid orange 
arrows are converged CCSD and A-CCSD(T) results from [381 . 

M. 



SRC and CC energies is, of course, a manifestation of how 
the lowering of the resolution scale improves not only the 
convergence w.r.t. the model space size, but also the con- 
vergence of the many-body expansion itself, as observed 
in the discussion of Fig. 2] (also see (l3j). 

Another interesting observation is that the IM-SRG(2) 
and A-CCSD(T) energies are surprisingly close, as pre- 
viously observed in calculations with SRG-evolved NN 
Hamiltonians in [f|. As discussed in Sect. Ill Dl the IM- 
SRG can be understood as a re-summation of the pertur- 
bation series for the ground-state energy, and the same is 
true for CC. We have argued that the IM-SRG(2) energy 
contains (at least) the complete third-order expressions, 
similar to CCSD, while A-CCSD(T) is complete through 
fourth order Q. Based on this assessment, and a su- 
perficial similarity between the CC cluster operator and 
the IM-SRG generator 77, one would have expected IM- 
SRG(2) and CCSD to give similar results. There are, 
however, differences between the two approaches: for in- 
stance, the effective Hamiltonian is always Hcrmitian in 
the IM-SRG, which implies that the method has more in 
common with unitar y C C (UCC) approaches rather than 
the traditional CC [42|. UCC converges more rapidly 
to exact Configuration Interaction results than its tradi- 
tional CC counterpart at a given truncation level because 
it contains additional many-body contributions, see, e.g., 
H|. The proximity of the IM-SRG (2) and A-CCSD(T) 
results in our calculations suggests that the same is true 
for the IM-SRG. 

For the NN + 3iV-induced Hamiltonian, we can ex- 
tend our systematics to A — > 00, i.e., the initial point of 
the free-space SRG evolution, by performing calculations 
using the bare NN Hamiltonian, with no induced 3N 
interactions. The results of such calculations arc shown 
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in Fig. [HI along with essentially converged CCSD and A- 
CC SD( T) results in CM ax > 18 spaces by Hagen et al. 
[H, Ull. The difference between these results and our 
eiviax = 14 CC results illustrate the slower convergence 
w.r.t. eMax for the bare NN Hamiltonian in comparison 
to the SRG-evolved Hamiltonians (cf. Sect. IIII( ). 

The calculations with the bare NN Hamiltonian show 
that the IM-SRG(2) energy lies below the exact NSCM 
result for 4 Hc, which is possible because the method 
is not variational. In fact, the IM-SRG(2) energy lies 
below the A-CCSD(T) energy, i.e., the deviation from 
the exact result is larger than for A-CCSD(T). The IM- 
SRG(2) energies for 16 0, and 40 Ca also lie below their A- 
CCSD(T) counterparts. There are different mechanisms 
which might explain this observation: The IM-SRG(2) 
might contain diagrams beyond A-CCSD(T) that lead to 
an overestimation of the binding energy, or IM-SRG(2) 
might lack diagrams compared to A-CCSD(T) that can- 
cel other contributions. A detailed diagrammatic anal- 
ysis of the IM-SRG which should clarify this issue is in 
progress, and will be a subject of a forthcoming publica- 
tion m. 

It might be tempting to consider the differences be- 
tween IM-SRG and CC energies for 4 He and the exact 
NCSM result for the bare NN Hamiltonian as a measure 
of uncertainty for the different many-body methods, be- 
cause it only depends on eMax and the many-body trun- 
cation. We emphasize that this is not possible, because 
the quality of the many-body truncation manifestly de- 
pends on the resolution scale A, as discussed in the con- 
text of Fig. [5] above. While the uncertainty due to the 
use of truncated methods like IM-SRG and CC is reduced 
for soft interactions, the omission of induced AN, . . . , AN 
interactions from the free-space SRC evolution, the omis- 
sion of the normal-ordered 3iV interaction, and the Esm&x 
cut in the present work all introduce additional uncer- 
tainties, which are hard to disentangle (also see (l3l|). 

Let us focus on the NN + 3A^-induccd Hamiltonian 
first. The IT-NCSM ground-state energies for 4 He and 
16 shown in the left column of Fig. [5] exhibit almost no 
A-dependence in the range from A = 1.88 to 2.24 fm _1 , 
which indicates that omitted induced AN, . . „ AN inter- 
actions are irrelevant for these light nuclei [j, [28| . The 
N02B approximation, i.e., the omission of the normal- 
ordered 3N term in the initial Hamiltonian (see Sect.HI|. 
causes an overestimation of the CCSD binding energ y by 
2% in 4 He, and 1% in the other nuclei studied here (l3j |. 
This overestimation is independent of A. If we assume 
that this remains true for IM-SRG(2) and A-CCSD(T) as 
well, the variation of the IM-SRG and CC energies with A 
must be caused by a combination of the many-body trun- 
cation and the £3 Max cut. Both of these effects become 
weaker as A decreases: The energy gain from including 
additional many-body terms is reduced as the conver- 
gence of the many-body expansion is improved, and the 
artificial energy gain that is caused by the omission of 
repulsive 3A~ matrix elements beyond the -E^Max cut is 
reduced as well. 



In 4 He, the deviation between the IM-SRG(2) and IT- 
NCSM ranges from 2.5% at A = 2.24 fin -1 to 1% at 
A = 1.88 mi" 1 . The effect of the £ 3Ma x cut is negligi- 
ble here, causing changes in the energy in the single- keV 
range. Thus, the deviation is dominated by the 2% con- 
tribution of the N02B approximation, which amounts to 
an overestimation of the exact NCSM energy by about 
500 keV. If we take this into account by shifting the IM- 
SRG and CC results upward, the IM-SRG(2) energies lie 
within a band of 350 keV around the NCSM result, which 
leads us to estimate the uncertainty associated with the 
truncation to be on the order of 1-1.5%. 

Size-extensivity of the IM-SRG(2) implies that the rel- 
ative truncation error, defined as the deviation from the 
exact NCSM result, should be of similar size in heav- 
ier nuclei. Adding the uncertainties associated with the 
N02B approximation and the i?3Max cut, which are both 
about 1% (see Sect. IIII|) . the total estimated uncertainty 
for 16 and 40 Ca is about 3.5%. For 16 0, this estimate 
agrees very well with the largest deviation between IT- 
NCSM and IM-SRG(2), which is 3.6%. We emphasize, 
however, that these estimates rely on the assumption that 
the effect of the N02B approximation is the same for IM- 
SRG(2) and A-CCSD(T) as for CCSD. 

For the NN + 3iV-full Hamiltonian (right column of 
Fig. [5|), we face the additional complication that the 
IT-NCSM energies for 4 He and 16 O are themselves A- 
dependent, which indicates that induced AN, . . . , AN in- 
teractions are no longer irrelevant. The underestimation 
of the experimental 4 Hc binding energy in the IT-NCSM 
calculations with the 7ViV+3./V-full Hamiltonian is caused 
by the omitted induced AN interaction (the fit of the 
low-energy constants was performed with the bare 3N 
interaction). Since the NN part of the Hamiltonian is 
unchanged from the NN + 3N- induced case, we can con- 
clude that AN, . . . , AN interactions must be generated by 
the SRG evolution of the initial 3 A interaction [28| . 

In He and 16 O, the relative differences between the 
IT-NCSM results and the IM-SRG and CC energies are 
reduced in comparison to the NN + 3Af-induced calcu- 
lation, to about 1.5% and 3%, respectively. Interest- 
ingly, the IM-SRG(2) energies now exhibit a weaker A- 
dependence than in the NN + SA^-induced calculations. 
This is consistent with the A-CCSD(T) results (also see 
[HI) and can be explained by the influence of the -E^Max 
truncation and the omitted induced AN, . . . , AN interac- 
tions in the case of the NN + SA'-full Hamiltonian. In 
4 He, the omitted induced interactions are attractive and 
their effect on the IT-NCSM energies grows only weakly 
as A decreases. For 16 0, on the other hand, the reduction 
of the IT-NCSM ground-state energy with decreasing A 
indicates that the net effect of the omitted interactions 
is repulsive. As mentioned above, for IM-SRG and CC 
we are also omitting parts of the the repulsive 37V inter- 
action by introducing the -E^Max cut. The A-dependencc 
of the two effects is reverse: With decreasing A, the ef- 
fect of omitted AN, . . . , AN interactions increases, but 
the contribution of matrix elements beyond the E^Max 
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FIG. 7. IM-SRG(2) ground-state energy per nucleon of closed- 
shell nuclei for NN + 37V-induced (top) and NN + 3iV-full 
Hamiltonians (bottom) at different resolution scales A. Ener- 
gies are determined at optimal Kl for eMax = 14. Experimen- 
tal energies (black bars) are taken from [371 ] - 



cut decreases. This leads to an apparent reduction of the 
A-dependence compared to the NN + 37V-induced calcu- 
lations. 

We conclude the discussion of our results by giving 
an overview of our IM-SRG(2) ground-state energies in 
Fig. the corresponding numbers are listed in Table [U 
With the NN + 3A^-full Hamiltonian, we obtain a good 
agreement with experimental binding energies in the s- 
and p-shell nuclei, while those of the sd- and p/-shell 
nuclei are overestimated by 0.6 — 0.8 MeV per nucleon, 
keeping in mind the previously discussed uncertainties of 
3-4%. A striking feature of Fig. [7] is the improvement 
of the isotopic trends of the ground-state energies with 
the inclusion of the initial 3iV force, in particular for the 
calcium isotopes. Since the effect of a A- variation on the 
ground-state energies is mostly uniform, this improve- 
ment appears to be a robust feature. 

For the NN + 3iV-full Hamiltonian, we also find that 
the aforementioned cancellation which reduces the A- 
dependence of the energies (cf. Fig. [5]) extends to all nu- 
clei studied in this work. This stability over a wide mass 
range is quite contrary to the growing A-dependence with 
A which is obtained without the initial 3iV Hamiltonian. 
It will be interesting to see if these systematica extend to 
higher masses, but an extension of our calculations be- 
yond the p/-shcll will require significant developments in 
the practical handling of the 3N matrix elements. 

Table U also contains energies which are extrapolated 
to an infinite HO space, using the extrapolation method 
proposed in [29| (also see Fig. [3]). A big advantage of this 
method is its ability to perform a simultaneous extrapo- 
lation for all calculated energies with different e^iax arL d 
Hfl, and thereby better constrain the fit model than sepa- 
rate extrapolations at each Ml. We find the extrapolation 







NN + 3N induced 


NN + 


3N full 








E[MeV] 






A^m" 1 ] 


eMax = 14 


extrap. 


eMax = 14 


extrap. 


4 He 


1.88 


-25.65 


-25.65 


-28.32 


-28.33 




2.00 


-25.75 


-25.75 


-28.38 


-28.38 




2.24 


-25.98 


-25.97(1) 


-28.51 


-28.51 


16q 


1.88 


-122.08 


-122.02(7) 


-130.73 


-130.6(1) 




2.00 


-122.80 


-122.71(9) 


-130.57 


-130.5(1) 




2.24 


-124.30 


-124.1(2) 


-130.47 


-130.3(2) 


24q 


1.88 


-155.76 


-155(1) 


-169.49 


-169(1) 




2.00 


-157.15 


-156.8(9) 


-169.47 


-169(1) 




2.24 


-159.66 


-159.4(6) 


-169.36 


-169.1(9) 


40 Ca 


1.88 


-354.14 


-354.1(10) 


-376.63 


-376(1) 




2.00 


-358.50 


-358.5(8) 


-376.26 


-376.1(8) 




2.24 


-365.50 


-365.7(8) 


-374.87 


-374.9(7) 


48 Ca 


1.88 


-411.34 


-409(4) 


-453.55 


-451(4) 




2.00 


-418.23 


-417(4) 


-454.36 


-453(5) 




2.24 


-428.74 


-429(5) 


-454.17 


-453(5) 


48 Ni 


1.88 


-337.74 


-336(3) 


-377.72 


-375(4) 




O AA 
Z.UU 


-344.30 


-343(4) 


-378.36 


Q'7'7/' A \ 

—61 ly<i) 




2.24 


-354.26 


-354(5) 


-377.60 


-377(5) 


56 Ni 


1.88 


-449.35 


-453(3) 


-529.62 


-530(4) 




2.00 


-460.10 


-464(5) 


-532.79 


-534(4) 




2.24 


-476.61 


-481(7) 


-534.48 


-536(6) 



TABLE I. Ground-state energies for closed-shell nuclei from 
IM-SRG(2) calculations with NN + 37V-induced and NN + 
3iV-full Hamiltonians at various values of A. Ground-state en- 
ergies for eMax = 14 are taken at the optimal value of Kl. Ex- 
trapolated energies are obtained using the method proposed 
in [29J. The uncertainties given in parentheses are solely due 
to the extrapolation, and determined by the largest difference 
between the eMax = 14 data for all Kl and the extrapolated 
values. Extrapolation uncertainties in 4 He which affect only 
the third decimal place are not listed. 

works well for the NN + 3iV-induced and NN + 3iV-full 
Hamiltonians, and all studied nuclei. As an error esti- 
mate, we include the maximum value of the scatter of 
eMax = 14 results around the extrapolated values E x 
in Tab. HI Since we are limited to Ml > 24 MeV in the 
calcium isotopes as well as 48 Ni, the extrapolation uncer- 
tainties are somewhat larger than for the lighter nuclei. 
Although the Ml = 24 MeV data for the A = 48 nuclei 
exhibit signs of the artifact which is caused by the trun- 
cation of the free-space SRG model space used for the 3 AT 
evolution (see Sect. IIIII and also Fig. [3]), the exclusion of 
these points from the fit data set does not affect the ex- 
trapolation within its uncertainties. For 56 Ni, however, 
the effect is more pronounced, hence we only include data 
points for Ml > 28 MeV when we determine the extrap- 
olated energies in this ca.se. 

V. CONCLUSIONS 

We have used the IM-SRG to perform a systematic 
study of closed-shell nuclei using chiral NN + 3N Hamil- 
tonians in the normal-ordered two-body approximation. 
Our calculations can be carried out in large HO bases 
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of up to 15 major shells due to the polynomial scaling 
of the IM-SRG with the single-particle basis size. We 
have demonstrated that the SRG evolution of the initial 
NN or NN + 3N Hamiltonian to a lower resolution scale 
is under control: Uncertainties caused by the omission 
of induced AN, . . . , AN forces are outweighed by the im- 
proved convergence behavior of the evolved Hamiltonian 
in many-body calculations, which reduces truncation and 
extrapolation uncertainties. 

With an initial Hamiltonian consisting of the N 3 LO 
NN interaction with cutoff A = 500 MeV, and an N 2 LO 
3iV interaction with a lowered cutoff A = 400 MeV, we 
achieve a very good agreement of the IM-SRG(2) ground- 
state energies with experimental data in light nuclei, 
while calcium and nickel isotopes are over-bound by 0.6- 
0.8 MeV per nuclcon, with minimal dependence on the 
SRG parameter A. Our IM-SRG ground-state energies 
are also in good agreement with results from quasi-exact 
IT-NCSVL as well as CCSD and A-CCSD(T) calculations 
HH EH HI] ■ The similarity of our results to the latter 
raises interesting questions regarding the diagrammatic 
content of the IM-SRG method when it is interpreted 
as a non-perturbative re-summation of the MBPT series. 
This aspect will be explored in depth in a forthcoming 
publication plj . 

Improving on the IM-SRG(2) truncation introduced in 
this work is straightforward: As a next step, we plan 
to implement the IM-SRG(3) flow equations, which con- 
tain explicit three-body operators. While the computa- 
tional treatment of 3iV interactions is expensive and re- 
quires additional approximations like the E^Max cut we 
discussed in the present work, it should be feasible to 
perform calculations in light nuclei to establish the con- 
vergence behavior of the IM-SRG hierarchy of trunca- 
tions. Effects in heavier nuclei can then be inferred using 
size-extensivity arguments. 

In the near term, we can use the existing technology 



to study NN and 3N Hamiltonian at consistent cutoffs 
and orders of the chiral expansion, in particular the re- 
cently completed N 3 LO interactions 0, |45[. We also 
plan spectroscopic applications by using the final IM- 
SRG single-particle energies and residual two-body in- 
teraction as input for traditional Shell Model calculations 
0. In parallel, we will pursue the generalization and di- 
rect application of the IM-SRG to open-shell nuclei via 
a multi-reference formalism. There are indications from 
quantum chemistry that a unitary method like the IM- 
SRG is more robust than traditional CC approaches in 
such cases [42l.l46l.l47i. 
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